Gauge field induced by ripples in graphene. 
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We study the effects of quenched height fluctuations (ripples) in graphene on the density of states 
(DOS). We show that at strong ripple disorder a divergence in the DOS can lead to an ordered 
ground state. We also discuss the formation of dislocations in corrugated systems, buckling effects 
in suspended samples, and the changes in the Landau levels due to the interplay between a real 
magnetic field and the gauge potential induced by ripples. 



I. INTRODUCTION. 

The recent characterization of graphene sheets made 
up of a single layer of carbon atomsi^ has caused great 
interest. Their unusual electronic band structure, and 
the possibility of tuning the number of electrons lead to a 
number of interesting features, both from a fundamental 
perspective and because of its potential applications^ 4 -. 

The low energy electronic states of graphene are well 
described, in the continuum limit, by two decoupled two 
dimensional Dirac equations. The kinetic energy depends 
linearly on the lattice momentum. The perturbations 
due to some types of disorder, like topological lattice 
defects^, strains^, and curvature^ enter as an effec- 
tive gauge field. Curvature, strains and topological lat- 
tice defects are expected to exist in graphene, as experi- 
ments show a significant corrugation both in suspended 
samples 10 , in samples deposited on a substrat o 11 ' 12 , and 
also in samples grown on metallic surfaces 1 ^. 

The statistical properties of the two dimensional Dirac 
equation in a random gauge field have been extensively 
studie d 14 ' 15 ] 16 i 17 i 18 , in relation with the Integer Quan- 
tum Hall effect. It has been shown that the density of 
states develops a peak at zero energy when the disorder 
strength exceeds a certain threshold. Furthermore, be- 
yond a second threshold, there is a transition to a glassy 
phased where the local density of states is dominated by 
rare regions.—. 

In the following, we will apply the analysis in 17 to the 
specific case of graphene, where there are two Dirac equa- 
tions coupled to the same random gauge field. The model 
will be detailed in the next Section. We analyze in the 
following Section the statistical properties of the gauge 
field. The main results for the density of states are pre- 
sented in Section IV. Given a divergent density of states 
at the energy of the Dirac point, we consider the instabil- 
ities which may be induced by interactions. Alternative 
approaches to the interplay between gauge fields and in- 
teractions are given i n 20 ' 21 ! 22 , although they did not con- 
sider diverging densities of states. Sections VI analyze 
the related problem of the structural changes which can 
be induced by the same random strains which give rise to 
the gauge field, following the analysis in 2 ^. Section VII 
discusses a buckling transition in suspended graphene. 
Section VIII estimates the effects of ripples on density 



fluctuations in the quantum Hall regime and compares 
with recent data 2 ^. The main results of the paper are 
summarized in Section VI. 



II. THE MODEL. 

We analyze the gauge field induced by the height fluc- 
tuations of a graphene layer on a rough substrate. In such 
case one expects that the shape of the graphene layer is 
determined by a competition between the interaction of 
the layer with the rough substrate, which tends to im- 
pose a preferred height, and the elastic properties of the 
layer. A simple Hamiltonian which models these effects 
is: 
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where h(r) is the height of the graphene layer, ho(r) is the 
preferred height which can be assumed to follow closely 
the substrate height, and ^(r) are the two in- 

equivalent Dirac (iso)-spinors which can be defined in the 
graphene lattice, ?%-(r) is the strain tensor associated to 
the deformation of the graphene layer, given by: 
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The gauge vector acting on the electrons in eq.{T]) is re- 
lated to the strain tensor b y 25 ' 26 : 



A x (r) = - [u xx {r) - u yy (r)} 
A v(r) = -2-u xy {r) 



(3) 



where a ~ 1.4A is the length of the bond between neigh- 
boring carbon atoms, and (3 = C/3 where C is a constant 
of order unity and (3 = — d\og(t) / d\og(a) ~ 2 — 3 is 
a dimensionlcss parameter which characterizes the cou- 
pling between the Dirac electrons and lattice deforma- 
tions. Besides the symmetry arguments in 2 ^, we also as- 
sume that the coupling between the electrons and lattice 
deformations is through the modulation of the hopping 
between nearest neighbor it orbitals, t ss 3e V 27 i 28 . The 
rest of the parameters which determine the hamiltonian 
in eq.([T]) are the electron Fermi velocity, v-p = 3ta/2, 
the bending rigidity, k ~ leV, the in-plane elastic con- 
stants, A, fj, ~ leV A~ 2 . To model the interaction be- 
tween the graphene layer and the substrate we use a sim- 
ple quadratic expansion around the height ho(f) which 
minimizes the energy in the absence of elastic and elec- 
tronic energy, parameterized by a coupling g. The value 
of this parameter is less understood. Estimates based 
on the analysis of the electrostatic potential between 
graphene and SiO^ suggest that g ~ 1CP 2 — 10~ 1 meV 
A~ 4 . By comparing g and k, one finds that the pinning 
by the substrate dominates for length scales greater than 
l p ~ (n/g) 1 / 4 ~ 10A. The coupling g being strongly rele- 
vant, for I 3> l p it can be considered as effectively infinite 
and the graphene layer rigidly pinned to the substrate, 
h(r) « ho(r) for I 3> l p . Note that we assume here that 
effect of direct pinning of the in plane modes by the sub- 
strate are small and can be neglected. 



III. EFFECTIVE GAUGE FIELD. 

Expcriment o 12 i 29 suggest that the height of the 
graphene layer shows fluctuations of order h ~ 10A over 
scales I ~ 100A. Similar fluctuations have been observed 
in suspended graphene sheets^. We will assume that the 
effects of the height fluctuations can be described statis- 
tically over distances larger than Iq ~ 100A. We will then 
relate the correlations of the effective random gauge field 
to the (four point) correlations of the (random) height 
profile, h(r). The calculation is valid whether this pro- 
file arises from interaction with a static rough substrate 
(in which case for I ^> l p it directly relates to substrate 
correlations) or from any other mechanism such as in 
suspended graphene. 

An estimate of the magnitude of the effective random 
gauge field can be obtained by noting that the height 
change between neighboring lattice points is ~ aVh 
hence the distance change is ~ a(Vh) 2 and the mod- 
ulation in t is St ~ [3t(Vh) 2 . Hence the modulation 
in A is ~ St/vp ~ (3(S7h) 2 /a which yields an estimate 



for the variance of the random effective magnetic field 
B = [V x A) z : 
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where h ~ 10 A is the typical scale of the height fluc- 
tuations, as discussed earlier. Typical parameters allow 
for (7 = 0(1), within range of the transitions that we 
consider below. 

In order to perform a more detailed calculation of the 
effective gauge field acting on the electrons, we first com- 
pute the in plane displacement field, u(r) obtained by 
minimizing the elastic energy for a given realization of 
h(f), and then we estimate the strain tensor Uij(r). We 
define: 
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In terms of these quantities, the procedure described 
above gives for the effective magnetic field acting on the 
electrons: 

(3fc 2 -fc 2 )(A + Ai) 
B{k) = (A + 2^ X 

kyfxX\k) ~\- k x fyy{}z) 2k x kyf X y{lz) (7) 

We assume that the average properties of the height mod- 
ulations are described by translationally invariant corre- 
lation functions, in Fourier: 



(fij{Q)fki(q)) = Fijki(q) 



(8) 



and are of short range character, i.e. with a finite limit 
for ql < 1: 

Fijkl(q)\q->o = fSijSki + /' (SikSji + SuSji) (9) 

a tensor compatible with the hexagonal symmetry of the 
lattice parameterized by two dimensionlcss constants / 
and /'. Using eqs.([7]), ([5]), and we find for the corre- 
lations ((4]) of the effective magnetic field at small q: 
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where q x + iq y = qe %B , where the angle 8 is measured 
from a given lattice axis. The angular dependence of the 
correlation is consistent with the lattice symmetry; as we 
show below, only its angular average is relevant for the 
transitions. 

In Eq. ([5]) we have assumed that the 2 point func- 
tion of (V/i) 2 field has a finite q = limit. The exact 
bound §;\{dih(f)djh(f)}\ < \^n(r-r f )\ 1 / 2 implies that 
the roughness h ~ of the graphene sheet (hence of the 
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substrate if adsorbed) can be at most £ < 1/2 in the gen- 
eral case for Eq. ([5]) to hold. In a model with Gaussian 
distributed h the condition is £ < 1/4 and higher rough- 
ness would result in long range (LR) correlations in the 
disorder. Such LR correlations would presumably arise 
when quenching thermal fluctuations of a freely fluctuat- 
ing membrane (which has ( = 0.59^2.) although a precise 
estimate then requires taking into account non gaussian 
fluctuations, a non trivial calculation. Here we restrict 
to SR disorder and substrates such that Eq. (7) holds. 
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IV. ELECTRONIC DENSITY OF STATES. 

We analyze the electronic density of states near the 
Dirac point, E = 0, using the techniques discussed inil. 
The main difference with the cases considered there is 
the existence of two Dirac equations coupled to the same 
gauge field, with couplings of equal absolute value but 
opposite sign, see Ti e i ec in eq.Jl}. 

The bosonized version of the problem also contains two 
fields, which become two sets of coupled fields when the 
replica trick is used to integrate over the disorder. Fi- 
nally, we make the same variational ansatz as ir>i£. The 
simplest observable is the total density of states (DOS), 
which is self averaging and is just twice the DOS of a sin- 
gle Dirac equation (single layer problem as defined ir>ii) 
and behaves as: 



p(E) ~ E 2 I Z ~ 



with: 
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where if is a parameter which describes the kinetic en- 
ergy of the field in the bosonized version of the model, 
and, for the non interacting case which corresponds to 
the hamiltonian in cq.([T]) takes the value K = 1. The 
parameter a determining the exponent in (| L 2[) is found 
to be given by the angle average of 
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i.e. the strength of the random gauge field, consistent 
with the order of magnitude estimate (J5]). 

The change in the dependence of the exponent z on 
the strength of the gauge field, a, in eq.fTS]) is associ- 
ated with a phase transition in the disordered bosonic 
model. For a > a c = 2/K 2 the local DOS (averaged 
over regions of size up to L ~ |i?| -1 / 2 ) exhibits strong 
fluctuations and non gaussian tails (i.e its disorder av- 
erage being different from its typical value) due to the 
dominance of rare regions. Note that the divergence of 
the DOS at E = occurs at a = 1/K 2 < cr c , i.e. be- 
fore the freezing transition in the one layer problem as 
disorder is increased. 



FIG. 1: Sketch of the DOS in a finite size L region for z < 2. 
The thick line is a typical value, while the thin line represent 
the size of fluctuations that are enhanced below the energy 
L~ z . For z > 2 the DOS increases at small E and its typical 
value saturates at L~ 2+z . For z > 3 (frozen regime) the 
fluctuations become so strong that the average DOS grows 
as L~ 2+z where z — 1 + a > z. Such finite size fluctuations 
should be observable in tunneling experiments. 



The effect on the DOS of an additional smooth random 
scalar potential with variance 5, corresponding to local 
fluctuations of the chemical potential, induced by e.g. 
the substrate, has been discussed irtil. It leads to: 
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and exhibits a transition at a' c = l/(2K ) = <t c /4. This 
leads to a finite and non zero DOS at zero energy: 



p{E) - 5 {2 ~ z ^ z ' 



(16) 



a behavior which thus exhibits two distinct freezing tran- 
sitions. The divergence of the DOS at a = 1 (for K = 1) 
is in between these transitions. 

Although we will not study this aspect in detail here, it 
is also interesting to note that since the two Dirac equa- 
tions describing the two valleys (the two Fermi points) 
feel opposite random gauge fields, mutual correlations 
of the local DOS in the two valleys as measured by 
{pi(E, r)p2{E, r)) c are strong. They are found to exhibit 
a transition at a different value of disorder a = 1/(2K 2 ) 
as can be seen by a study analogous to the two layer 
model of Section IV B ofii 



V. INTERACTION EFFECTS AND 
ELECTRONIC INSTABILITIES. 

For sufficiently large disorder, a > 1, the density of 
states, and, as a consequence, the electronic compress- 
ibility, diverges at E = 0. The electron-electron interac- 
tion, or the interaction of the electrons with other degrees 
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FIG. 2: (Color online). Critical temperature as function of 
chemical potential. The parameters used are Wo = 200meV, 
lo = 10a, a = 1.4(2 = 2.4), and U = leV. The value of 
a implies an average height fluctuation h m 3.4A. The blue 
diamonds give the critical temperature when the transition 
is discontinuous. The green triangles are the values of the 
gap A, in Kelvin, at the transition temperature, in the region 
where A jumps discontinuously from zero to a finite value. 
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FIG. 3: (Color online). Approximate phase diagram, as func- 
tion of electron density and temperature, obtained with the 
same parameters used in Fig. [2] . The existence of a first order 
transition leads to a region where electronic phase separation 
is induced. 



of freedom, will lead to instabilities, which suppress the 
compressibility. 

Within mean field theory, the effects of interactions 
on the electronic band structure can be described as an 
external potential which must be calculated self consis- 
tently. A simple such potential which opens a gap at 
E = is the shift of the energy on one sublattice of 
the honeycomb structure with respect to the other. This 
shift can be associated to a spin or to a charge density 
wave, or it can be induced by phonons 3 ^ or by short range 
electron-electron interaction s 31 ' 32 ' 33 ' 34 . In the continuum 



model described here, it enters as a mass term: 
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where TL e i ec is defined in eq.([T]). The total electronic 
hamiltonian satisfies: 



(Helec + T~La) = {T~te 
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(18) 



where X is the four dimensional unit matrix, indepen- 
dent of spatial position, which acts on the space spanned 
by the four component electronic (iso)-spinors. The 

eigenvalues of the electronic hamiltonian satisfy e l ° t2 = 
e„ + A 2 , where e n is an eigenvalue of 7i e iec- As a result, 
the density of states associated to H'°* c , Pa(E), satisfies: 
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and, using the expression in eq. (|ll|) . we find: 
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where the energy scale Wq — v-p/lo is inserted so that for 
E Wq ^> A the value of pa{E) crosses over into the 
density of states of the clean system, p(E) ~ l^l/w 2 ,. 

The self consistent value of A is determined by the 
competition between the cost in energy associated to the 
formation of the gap, and the decrease in electronic en- 
ergy due to the reduction in the density of states at the 
Fermi level. Near the transition, A is small compared 
to the other energy scales of the model, and the energy 
required to create the charge or spin density wave can 
be expanded as function of A. For simplicity, we assume 
that the ordered phase is a spin density wave induced by 
the on-site Hubbard repulsion, U which breaks the sub- 
lattice symmetry and produces a gap A = U S/2 where S 
is the resulting polarization per site. The total energy is 
the sum of the kinetic energy and the gain in interaction 
energy obtained by inducing the polarization: 
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where 4 allows for spin and valley degeneracy and we 
allow for a possibly non zero Fermi energy Ep- The in- 
duced gap is given by minimizing the total free energy 
and with a change of integration variable: 
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We consider first the case Ep 
E el = -Aa*S Wo 



T = where F el -> 
EpA(E)dE. The integrand can be ex- 



panded for E A, where it goes as A 2 E^ 2 / Z ^ 2 . As a 
result, we obtain a contribution to the electronic energy 
^iE e i ec (A) ~ — (a/^ ) 2 A 2 /W . There is also a contribu- 
tion from the region E ~ A. This term in 2£ e / ec (A) can 

be written as <5 2 £ e i ec (A) - -(a/Z ) 2 A( 2 /^ +1 W ~ 2/z . 

The relative strength of the two terms discussed above 
leads to the existence of three regimes: i) 2/z — 1 > 
0. The electronic energy is determined by 5iE e i ec (A). 
Both the magnetic and electronic energy go as ~ 
A 2 , and, for U <C v-p/a the minimum energy is at 
A = 0. ii) For 2/z - 1 = 0, we find 8iE elec ~ 
-2(a/; ) 2 A 2 /Wolog(Wo/A). The magnetic energy is 
greater by a logarithmic factor, and there is an ordered 
phase, with A - w e- { - Wol oW Ual \ The problem be- 
comes equivalent to the Peierls analysis of the instability 
of a one dimensional metal, iii) For 2/z — 1 < 0, the lead- 
ing contribution is ^-E^iec- There is a magnetic phase 
with a gap: 
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We now analyze the way in which the magnetic phase 
which always exists for 2/z — 1 < is modified when 
Ep,T 7^ 0. In particular the order of the transition is 
determined by the sign of the 04 coefficient in the free 
energy expansion F = a 2 A 2 + 04 A 4 . Taking a 9a 2 |o 
on the right hand side of Eq. (|22|) yields 04, hence the 
simultaneous conditions a 2 = and 04 = determines 
a critical Ep, with Ep = a(z)T c whenever T c -C Wo, 
such that the transition changes from second order for 
Ep < E F , to first order in the region Ep > Ep where we 
have 04 < 0. We find numerically that cosh a{z) varies 
between 3.4 at z = 2 and 2 at z — > 00. 

A typical phase diagram with z > 2 is shown in Fig. [2] , 
where the value of the gap A c at the critical tempera- 
ture, in the region where the transition is first order is 
also shown. When the line of first order transitions is 
crossed, the electron density jumps discontinuously. For 
sufficiently large values of Ep we find that A c (Ep) > Ep. 
When the transition line is crossed in this region, the 
electron density in the ordered phase is zero. The phase 
diagram as function of temperature and electron density 
is shown in Fig. [3] . 



VI. FORMATION OF LATTICE DEFECTS. 
A. Unbinding of dilocations. 

As discussed earlier, ripples, e.g. due pinning to a 
rough substrate, induce in plane strains. If these strains 
are sufficiently large, it will become favorable to relax 
them by creating lattice dislocations. It is convenient to 
view the out of plane deformations as inducing quenched 
random stresses coupling linearly to the in plane strain 



tensor Uij via an energy density Putin . One can 



then apply the result of Ref. [23j for the threshold beyond 
which random stresses generate dislocations. 

Using eq.fl]), the random stress tensor field which 
renormalizes the fugacity of dislocations is: 



>.iy 



A 


(dh 


2 


\~dx~ 


A 


'(dh 


2 


\dx 


dh dh 


^ dx dy 



dh 
dy 

dh 
dy 



dh 



dh 
dy 



(24) 



We assume that the correlations of this field are given 
by: 

{(Jij(q)(T k i(-q)) | 9 _> = [cr\5ijS k i + <x M (5 ik 5ij + SuSj k )] 

(25) 

where the parameters and o\ + 2(7^ measure the 
strength of random shear stresses and comprcssional 
stresses, respectively. 

In presence of random stresses an isolated dislocation 
in a region of size L feels a random potential whose min- 
ima grow typically ~ —hxL. The logarithmic elastic en- 
ergy cost of creating a dislocation can then be overcome, 
and thus dislocations will proliferate at T = 0, when: 



. _ A(A + 2^)g A1 + A t 2 (q A + 2cr A1 ) _ a 2 
° ~ M 2 (A + ^) 2 - CTc ~16tt 



(26) 



where we have neglected the effect of screening of the 
elastic coefficients by disorder, which have been shown 
to be small 4 ^). To the same accuracy this formula 
holds for all T < T m /2 where T m = K a 2 /(16ir) is the 
KTHNY melting temperature of a pure 2d crystal, with 
A'q = 4jtx(/U + A)/(/i + 2A), while the threshold decreases 
as a c (T) = 4<7 C ^(1 - at higher T. For a > <r c 
and at T = the scale L above which dislocation first 
appear can be estimated as inM- and corresponds to the 



total energy cost K °£~ (1 — ^/<t/ct c ) ln(L/i ) + E c becom- 
ing negative. We have taken into account the dislocation 
core energy E c = E° + \n{l /a) at scale Z (E° de- 
notes the bare core energy). Because of logarithms this 
scale can be large hence it can alternatively be viewed as 
defining an effective size dependent threshold a c (L). The 
dislocation density above this scale can be estimated by 
arguments similar to^ 2 -. 

The quantities a\ and can be written in terms of 
the correlations of the function fy, given in eqs.([8]) and 
©: 



[ M 2 (/ + 2/') + A(A + 2 M ) (/ + /')] 



a » = J 2 » 1 



(27) 



Inserting this result in cq . (|2"6")) . and assuming that 
j3, A//i ~ O(l), we find that dislocations will proliferate 
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when the height correlations arc such that h 2 /(l^a) > 1, 
which is the same combination of scales which determines 
the existence of a divergence in the electronic density of 
states. 



B. Buckling into the third dimension. 

An effect not taken into account above is that dis- 
locations may buckle in the third dimension to lower 
their energy. For a free membrane (in the absence of 
a substrate) this occurs for scales larger than the buck- 
ling radius Rb, and below that scale the membrane re- 
mains flat and Coulomb gas logarithmic scaling holds. 
In principle, for a free membrane in presence of internal 
in plane random stresses, if Rb is large enough (values 
such as Rb ~ 10 2 k/(Kqo) are quoted in Ref4^), i.e. if 
Rb > Iq ^> a, the above energy estimate setting L = Rb 
can be used to determine the disorder threshold at which 
buckled dislocations would occur. However, if one takes 
into account the pinning of the height held to the sub- 
strate, the energy calculation of Ref.— ) remains valid for 
scales smaller than l p , but must be reexamined for scales 
larger than l p , a problem left for future study. 



C. Gauge fields associated to dislocations. 

Note, finally, that dislocation cores act on the electrons 
outside the core as vortices^ 2 - of flux $ = e$o/3, where 
<f>o is the quantum unit of flux (= 2n in our units), and 
e = ±1. Hence, the existence of dislocations will increase 
the random field due to elastic strains considered so far. 
Given a set of dislocations at position f n and Burgers 
charges b n the resulting effective magnetic field can be 
written B{r) = ($o/3)n(r) where n{r) = ^2 n £ n S(f — 

f n ) and the signs are given by e = 2b ■ ciimod27r. If 
positions and signs were chosen uncorrelated (such as in 
a quench from infinite temperature) it would result in a 
LR correlated random gauge field, i.e Cs{q) ~ ^d^ 2 
at small q in Q, where d is the mean distance between 
defects 4 ^ 

This procedure however leads to Burgers charge fluc- 
tuations growing as ± ~ L in an area L 2 hence a very 
large elastic energy, L\nL. If the system can relax, this 
energy is screened and the result is a finite parameter 
a as defined in In cases where the dislocation den- 
sity is not very small it can be estimated from a Dcbyc- 
Hiickel theory. One non-equilibrium example is a quench 
of a pure crystal from a (moderate) temperature Tq > 
T m to low temperature in which case (n(q)n(— q)) = 
T Q q 2 /{E Q c q 2 + K a 2 ) hence a = T Q <t> 2 /(9nK a 2 ). Fur- 
ther relaxation of ~ In L energy would then occur. An- 
other example is the distribution of dislocations induced 
by the ripples as in (|26|) . with a > a c . Then one 
estimates 4 ^ (n(q)n(-q)) = jq 2 aK 2 a 2 / (E®q 2 + K a 2 ) 2 
hence a = <j^>q/ (36Tra 2 ). Very near the transition Dcbyc- 



Hiickel docs not apply as a vanishes at a c proportionally 
to the density of dislocations. 



VII. RIPPLES IN SUSPENDED GRAPHENE 

Finally we discuss a possible source for ripples in sus- 
pended graphen o 10 i 43 ' 44 . Upon etching a preexisting 
rough substrate, the rippled graphene sheet would tend 
to relax to a flat configuration with higher projected area. 
This however may be precluded if the sheet is pinned at 
its boundaries. Indeed it is known that fixed connectivity 
membranes exhibit a buckled state when constrained at 
their boundaries by a fixed frame of projected area Af 
smaller than the equilibrium area of the unconstrained 
membrane A. As discussed irt 4 ^ it results in an additional 
comprcssional energy term of the form r J d 2 r^2 i Uu and 
hence implies that the energy of flexural modes becomes, 
to lowest order \ J d 2 f[K(S7 2 h) 2 + r(V/i) 2 ]. In the buck- 
led phase, Af < A, r < 0, an instability thus develops 
at scales larger than ^ ~ (k/|t|) 1//2 . This phase can 
be described as a non homogeneous mixture of pure flat 
phases with different orientations. For arbitrary bound- 
ary conditions, it is expected to be non trivial since, con- 
trary to a one dimensional rod, in a polymerized mem- 
brane the in-plane modes cannot fully relax the flexu- 
ral constraints, i.e. the transverse part of the flexural 
strain tensor, P?(y)(dihdjh), cannot be relaxed by the 
in-plane strain field. While demonstrating the instability 
is simple, the full calculation of the resulting shape re- 
quires consideration of non linear terms and is difficult. 
The problem of relaxation from a randomly rippled con- 
figuration with a fixed frame constraint deserves further 
study, in particular the question of whether there is some 
memory of the initial ripple pattern. 

Note that one may consider, alternatively, uncon- 
strained boundary and apply a tension — / J d 2 rdiUi. 
The buckling transition 4 ^ has been mostly studied on the 
side where the sheet is stretched (the effective 77? — ► + ). 
It was found that (A f - A) /A ~ t r ~ |/| 1 / ,5 sign(/) and 
the correlation lengths ^ ~ £ u ~ l/l - "/ 5 for flexural and 
phonon modes at small /. While entropic effects produce 
non trivial values for these exponents, these may be ob- 
servable only at large scales, and at intermediate scales 
mean field values 5 = 1, v = 1/2 (discussed above) are 
appropriate. It would thus be interesting to study the 
other side of this transition. 



VIII. BROADENING OF LANDAU LEVELS IN 
PRESENCE OF A REAL MAGNETIC FIELD 

In presence of a real magnetic field B we expect that 
the Landau levels will be broadened by the effective field 
B rip due to the ripples. Alternatively, the local electronic 
density corresponding to N full Landau levels is fluctu- 
ating according to n(r) = [B + B r i p (r)]N/<f>o, where <j>o 
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is the flux quantum. Such density fluctuations were re- 
cently measured^ showing Sn = ±2.3 • 10 n cm~ 2 at a 
field of 11T for N = 2, 6, 10. In this section we estimate 
the contribution of the random gauge field due to ripples 
to these density fluctuations. 

Consider the density n measured on a length scale L, 
and its probability distribution, near an average density 



BN/ipQ. We assume first > Ib\N, where Ib = ^<po/B 
is the cyclotron radius and IbVN estimates the size of an 
orbit in the N-th Landau level. Each Landau orbit has 
then a random shift ±B r i p (r) where the ± corresponds 
to the K and K' valleys that feel opposite gauge fields. 
The density distribution is 

P(n)= f J2 S ( n - BN ^ 0±B "P^ N /M d2r / 2L2 

Jr<L ± 

(28) 

The average is (n) = BN/tpo while the variance is: 



N 2 

(Sn 2 ) = -2^4<[/ 



N 2 



r<L 



= 4^£4«f A{u)du] 2 ) 
(29) 

where B r i p = [V x A] 2 0o/27r and A is the random gauge 
field considered in the previous Sections (with the appro- 
priate change in units) and the contour encloses the area 
of measurement. To estimate the variance in (|2"9")l we use 
Eq. §5§ with a cutoff exp (— q 2 l 2 /2). In real space it cor- 
responds to (Ai(r)A,-(r)) - (cr/2Z 2 )cxp[-(r-r') 2 /2^] 
where we neglect*^ the transversality constraint on A. It 
yields: 



(Sn 2 



N 2 a 
47r/ L 3 



(30) 



The experimentally more relevant case is Iq < IbVN. 
In this case we argue that A(r) can be replaced by its 
average within a Landau state, so that an average ripple 
field is 



B av (v) 



1 



2ttN1 2 b 



d 2 r B rip (r )e 



-(r-r ) 2 /2JVil 



(31) 



iV 



and its Fourier transform is B av (q) 
B r ip(q) exp(— q 2 Nl 2 B /2). This replaces Iq — > Ib 
in Eq. (|30|) . and identifying L with the tip size lu p in 
the experiment 2 ^, we obtain, 



(Sn 2 



_/V 3 / 2 



(32) 



tip 



UsingSi l w lOOnm and l B « lOnm, Eq. (021) yields 
numbers consistent with the experiment, except for the N 
dependence. Note that an even weaker dependence in N 
(Sn ~ N 1 / 4 ) is obtained if one assumes that L = IbVN 
is the only averaging scale. 

It is also interesting to estimate the energy broad- 
ening <5ej\r of the Landau levels, which in the absence 
of ripples have energies ejy = vp\[2e\/ BN' with N = 
AN' + 2. The field associated with the ripples changes 



locally the energy of the Landau levels, which become 
e N = v fV2c\/ (B ± SB r i P )N', where SB rip is the aver- 
age value of B r i p in the region occupied by the Landau 
level. A similar calculation then yields the estimate for 
N > 2 



Se N > 



v F 



1/2 



AT-i/4 



(33) 



in the regime IbvN > Iq. 

Finally we note, that the N' = level has no broad- 
ening at all. This remarkable result is obtained by fac- 
torizing the N' = eigenstates of the free Dirac system 
in a magnetic field with the well known zero energy so- 
lutions of the random gauge proble m 14 ' 17 . This set has 
the proper Landau degeneracy and is therefore an exact 
solution for the zero energy Landau level with random 
gauge. 



IX. CONCLUSIONS. 

We have analyzed the effect of random gauge fields 
on the electronic structure of corrugated graphene. We 
find that the local density of states diverges at the Dirac 
energy E — 0, as p(E) cx E 2 ^' 1 , with z > 2, for suffi- 
ciently strong disorder. The scale of height fluctuations, 
h should satisfy f3h 2 /(la) > 1, where /3 ~ 1 — 2 gives the 
coupling between the electrons and the lattice strains, I 
is the typical spatial scale of the disorder, and a is the 
lattice constant. 

A divergence in the density of non interacting density 
of states implies the existence of instabilities in the pres- 
ence of electron-electron interactions. We have analyzed 
the possibility that a gap will open at low temperatures, 
depicting the low energy density of states. We have found 
a first order transition to an ordered state at large Ep. 
This discontinuous transition, in turn, implies electronic 
phase separation. 

When the strains which induce the gauge potential are 
sufficiently strong, they can lead to an instability, and the 
formation of lattice dislocations. This change takes place 
for C(X, ix)h 2 /(la) > 1, where C(A, /z) ~ 1 is a dimension- 
less parameter which depends on the elastic constants of 
the material. 

We have described the main features of the buckling 
instability which may arise in suspended systems under 
compression. Finally, we analyze the changes induced 
in the Landau levels induced by a magnetic field by the 
gauge potential associated to ripples and show correspon- 
dence with experimental data" 

Our analysis is consistent with previous work on the 
changes in the electronic density of states in graphene in 
the presence of ripples 3 ^ (see also^i). A transition to a 
state magnetically ordered in highly disordered systems 
agrees with the observation of magnetism in irradiated 
graphite samples^. The existence of charge inhomo- 
geneities, due to electronic phase separation, can help 
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to explain the observation of charge puddles when the 
Fermi energy is close to the Dirac energy^. 
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